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Searches for gravitational wave signals which do not have a precise model describing the shape of 
their waveforms are often performed using power detectors based on a quadratic form of the data. 
A new, optimal method of generalizing these power detectors so that they operate coherently over 
a network of interferometers is presented. Such a mode of operation is useful in obtaining better 
detection efficiencies, and better estimates of the position of the source of the gravitational wave 
signal. Numerical simulations based on a realistic, computationally efficient hierarchical implemen- 
tation of the method are used to characterize its efficiency, for detection and for position estimation. 
The method is shown to be more efficient at detecting signals than an incoherent approach based 
on coincidences between lists of events. It is also shown to be capable of locating the position of the 
, source. 
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I. INTRODUCTION 

Six kilometer-scale laser interferometers designed to observe gravitational waves (GW) with unprecedented sensi- 
^vq ■ tivities should complete or approach the end of their commissioningin the year 2003. Three are operated in North 
| America by the LIGO Laboratory pj, two in Europe by the Virgo and the GEO600 Q projects, on one in Asia 
by the TAMA300 t 4iJ project. A collaborative analysis of the data collected by these instruments provides the best 
OO ' prospects for detecting and analyzing GW events of astronomical origin. 

The focus of this article will be on "bursts" of gravitational radiation, loosely defined as GW of duration of the 
order of a few seconds at most, and present in a frequency range overlapping at least partially with the bandwidth of 
the interferometers (10 Hz < / < 1 kHz). Other types of signals that will not be discussed here include continuous 
^ GW from rotating neutron stars, and a stochastic background of GW of cosmological origin. 

Arguments based on the astrophysics and on the dynamics of the sources of GW bursts show that the detection 
of these signals will be challenging, as the signals will be deeply buried in the instrumental noise Consequently, 
^P" a significant research effort is currently on-going to develop and study efficient algorithms for the detection and the 
characterization of the elusive GW signals. An important fraction of the literature on the subject concerns signals 
with a precisely known form 6] . The knowledge of the signal allows the construction of a phase coherent filter (the 
Wiener or matched filter) which is known to be optimal for signal detection. Only the coalescence of compact binaries 
and possibly the ringdown of excited black holes should be detectable using matched filtering. For the particular case 
of compact binaries, it is known that a coherent analysis using data from all the interferometers of the international 
network will improve the detection prospects noticeab ly m |8j , although the computational cost of such an analysis 
might be prohibitive 0. In addition, it was shown in |ldj| that the use of the Advanced LIGO detectors and of the 
Virgo interferometer cooperatively might allow the localization of the GW source with enough accuracy to permit its 
observation with electromagnetic instruments, thus complementing with information about the thermodynamics of 
the source the information on its dynamics provided by the GW. 

The remainder of this article will be concerned only with GW signals that are not known with enough precision to 
allow matched filtering. The algorithms that have been proposed in the literature to detect these signals fall into two 
general categories: time-domain filters, and power detectors. Time-domain filters fri| rely on the development of a 
small bank of linear filters which are expected to cover relatively well the space of possible GW signals. They offer 
the advantages of speed, simplicity, and possibly ease of interpretation, but might lack the robustness and efficacy of 
power detectors. Only the latter will be discussed here. The power detectors threshold on some non-linear measure 
of the data, often constructed from a time-frequency representation of the signal 0, 0, 0, 0, 0, . They have 
been shown to be optimal for the detection of signals with especially poor waveform descriptions. 
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All power detectors were explicitly designed and implemented to process data from the different interferometers of 
the world-wide network independently. Under this mode of operation, it is expected that event lists are generated 
individually from the data stream provided by each interferometer, and are later compared to form coincidences 
based on temporal, frequency, or more general information. This incoherent approach should not yield the maximum 
efficacy, in part because GW bursts in individual interferometers have to be rather loud to register with the power 
detector and to give accurate estimates of their start time, duration, frequency band, amplitude, etc., all of which 
might be needed by the coincidence gate. The alternative is to combine all data streams first, and then run a burst 
detector on the synthetic data stream so produced. I implement this coherent approach as a generalization of the 
power detectors developed to date for single interferometers, by calculating the optimal way to combine any number 
of interferometer data streams into a single time-series, such that when this time-series is fed to a single interferometer 
power detector, a larger signal-to-noise ratio is obtained than for any other combination of the data streams. This 
brings to the already implemented and well-characterized power detectors the benefits of a network coherent analysis, 
which include improved sensitivity, and the ability to precisely locate the source position on the sky. 

A. Summary of Results 

The coherent power filter (CPF) algorithm presented in this paper involves the following steps: 

1. A point on a grid in parameter space is chosen, where the parameters are the source angular position, and two 
numbers describing the plus and cross polarizations of the GW signal. One of these numbers is the ratio of the 
power in the cross polarization to the power in the plus polarization, and the other is the amount of overlap 
between the two polarizations, approximately measured by integrating the product of the two waveforms over 
time. 

2. Given the source position and the network topology, the data streams from all interferometers are time-shifted 
to align the GW signals to a common origin in time. 

3. Every data stream is multiplied by a scale factor, and all the data streams are added together to form a synthetic 
time-series. The choice of the scale factors depends on the network topology and on the four parameters chosen 
in Step 1. 

4. The synthetic time-series is processed by a power detector, and the power measurement is recorded. 

5. If all points of the grid in parameter space have been visited, the algorithm exits. If the maximum of the power 
measurements exceeds a certain threshold, a detection is announced, and the parameters that gave the largest 
power measurement are returned as an estimate of the source parameters. 

6. Back to Step 1. 

A more detailed discussion of this algorithm is presented in section IIIII 

The scale factors in Step 3 are chosen so that the value of the signal-to-noise ratio (SNR) is maximized. It is quite 
significant that only two parameters in addition to the source position are required to perform this maximization. 
Geometrically, this can be understood by realizing that the signals in all interferometers, after being properly time- 
shifted, are linear combinations of the two polarization waveforms, and therefore lie in a hyperplanc spanned by these 
two polarization waveforms. Consequently, a knowledge of the ratio of the lengths of the two polarization waveforms 
and of their angle with respect to each other, together with a knowledge of the beam-patterns of all interferometers, 
is sufficient to determine the signals in all interferometers, up to an overall scale factor, and up to the orientation of 
the hyperplane. However, neither of these two pieces of information are needed to calculate how to linearly combine 
the signals to get maximum power in the synthetic time-series. 

The validity of the geometrical picture and the actual conclusion that only four parameters are required to perform 
the SNR maximization depends critically on the right choice of the position dependent time-shifts. As shown in 
section IIIII in the specific case where the cross-correlation function of the plus and the cross polarizations does not 
have an extremum at zero time lag, there are no formal guarantees that the CPF algorithm will converge to the 
right source parameters. Physically, this results from possible interactions (or cross-terms) between the plus and the 
cross polarizations which cannot be properly handled by the coherent algorithm. This does not affect the detection 
performances of CPF, but in some cases is significant for source position estimations. As discussed in section IIIII a 
number of canonical sources do not satisfy this condition exactly, so that a careful study of the position systematic 
errors is needed. For the difficult case where the two polarizations are long monochromatic signals with a phase 
difference of a quarter of a cycle, it is shown in section IIII Al that a correction for this systematic error can be 
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implemented such that for ~ 25% of the sky the systematic error is negligible, while that for about 50% of the sky it 
is too large to allow any position estimation. This correction procedure only requires the additional knowledge of a 
quantity which is closely related to the characteristic frequency of the signal. 

The performances of the CPF algorithm are explored empirically in section lTvl through numerical simulations. All the 
experiments are limited to the three interferometer network (the HLV network) consisting of the LIGO interferometer 
near Hanford, Washington, the LIGO interferometer near Livingston, Louisiana, and the Virgo interferometer in Italy. 
The signal is short (1/16 s) and narrow band (25 Hz), and is assumed to originate from a position along the northern 
hemisphere normal to the HLV plane. All experiments are performed with the TFCLUSTERS 16] algorithm as the single 
interferometer power detector. A realistic, computationally efficient hierarchical implementation of the CPF algorithm 
is shown to offer better detection performances than a incoherent approach which uses only coincidences between 
events generated by independent TFCLUSTERS operating on the three interferometers. It is also shown that the CPF 
algorithm can be used to estimate the position of the source of GW. When the GW signal has four times more power 
in its plus polarization than in its cross polarization, roughly one quarter of all trials lead to a position estimate that 
is within one degree from the true source position, for signals with reasonable amplitudes. The ability to pinpoint the 
source location is debilitated by the misalignment of Virgo with respect to the LIGO detectors, the reduction of the 
signal-to-noise ratio, and the reduction of the ratio of the power in the plus and in the cross polarizations. 

II. NOTATION 

Let bold characters denote time-series; whether these time-series are continuous in time or discretely sampled will be 
immaterial in the following discussion. It is assumed that a GW is observed with a network of N independent detectors. 
Calibrated data corresponding to measurements of the GW strain in all detectors are denoted yi, i = 1,2, N. The 
noises rii are assumed to be additive, so 

y i = F % +(d,cl>,iP)T[A i (9,<l ) )]s + +F?"(6,<f>,^)T[A i (6,ct>)}s x +n t , (1) 

where s + and s x are the two polarizations of the GW signal, F* are the beam- pattern functions 5] of the i th 
detector, and T(A) denotes the time-shift operator; for time-series with continuous time, for instance, T(A)x(t) = 
x(t — A). The beam-pattern functions depend on the two angles describing the source position (the right ascension 
and the declination, denoted 8 and 4> respectively), and on the polarization angle if). The time-shift at the I th detector, 
denoted Aj, is the same for the two polarizations, and depends only on the source position on the sky. The frame in 
which s + and s x arc defined is irrelevant since the waveforms are not assumed to be known a priori; a rotation of 
that frame is equivalent to a change in s + , s x and tp. The parameters 0, 0, tp, s + , and s x , and the derived quantities 
F^~(8, <j>, -0), F* (8, </>, tp), and Ai(6, 0), will be used below to describe the parameters of a real source which is assumed 
to be present in the data, and which we are trying to detect. 

The scalar product between two time-series is denoted x ■ y. For time-series with continuous time, it is defined as 

/oo />oo 
/ x(t x )Q(t x ,ty)y(t y )dt x dty, (2) 
- oc J — oo 

and similarly for time-series with discrete time. The kernel Q can be viewed as a filter applied to the time-series in 
order to detect more efficiently a particular signal or to modify the character of the noise, for instance. The square of 
the norm of a time-series, also called its power, is denoted \x\ 2 , and is defined by \x\ 2 = x ■ x. 

The noise in each of the N interferometers is only assumed to be wide-sense stationary 0], i- e -> it does not have 
to be Gaussian or white. The noises can always be made zero mean and independent by linear filtering |19|. so 
E[yi ■ T(A)yj] = Ri(A) if i = j and is zero otherwise, for i?i(A) the autocorrelation of the noise, and E[-} denoting 
the expectation value of its argument. As usual, the Fourier transform of the noise autocorrelation function is the 
noise power spectral density. 

III. ALGORITHM 

The synthetic response of the network is denoted Y , and is a simple linear combination of the time-shifted individual 
detector responses: 

N 

Y = Y,Oi T (8i)Vi, (3) 

i=l 



4 



for some set of real coefficients a, and time shifts Si, i = 1,...,N, which are arranged in two vectors, ~ct and 5 , 
respectively. Note that the Si are the trial time delays used in the data analysis, and the algorithm defined below is 
used to estimate these delays so that they are close from the real delays Aj(0, 4>) corresponding to a source located at 
position (9, </)). The network power P is the estimate of the power in the GW signal, according to our norm definition, 
i.e., P=\Y\ 2 . 

The motivation behind this particular design is that it is the simplest generalization of the numerous power detec- 
tors for single interferometers described in the literature, which in many cases have already been implemented and 
characterized. In practical terms, a software code can be designed to compute the synthetic response for a network of 
interferometers, and these data can be fed to a power detector, as if they were data from a single interferometer, in 
order to measure the network power. The kernel of the dot product used for the computation of the network power is 
then determined by the single interferometer power detector used to process the synthetic response. Different power 
detectors are efficient for detecting different types of signals, so this generality of the synthetic response approach 
is very economical in terms of code development. Some single interferometer power detectors provide a non- linear 
measure of the power; this is not a serious limitation given the algorithm structure defined above, for the power 
measurements are all very nearly linear for detectable signals. 

Let Sij = S{ — Sj + Aj — Aj denotes the error on the estimated time-of- flight between detectors i and j. The network 
power can be expanded as 



P = C + v, (4) 



where the signal term is given by 

N 



C = J2 a i a j [F+F+R ++ (S ij ) + F+F*R +x (S l3 ) + 

F* F+R x + (Sij ) + F* F* R x x (Sij )] , (5) 

and where the noise term is given by 

N 

V=Y1 wATtfi + A i)(K +s + + p t Xs x) ■ t ( s j + A j) u j + T ( S * + A *K • T(6j + Aj-Jnj]. (6) 

The signal correlation functions are given by 

Rijik - tj) = T(t i )s l ■ T(tj)sj, (7) 

for i, j = + or x . 

The signal-to-noise ratio p is defined as 



where the expectation of the noise can be rewritten as 

N 

E[ V ] = J2*>1 (9) 

i=l 

where the noise variance is of = Ri(0). For fixed noises and signals, the signal-to- noise ratio depends only on 7? and 
S . It can therefore be maximized for a given choice of the source parameters (9, <j), ip, s+, s x ) by varying ~ct and S . 
Let ~rf m (8\ <f)\ ip', s' + , s' x ) and S m (6' , <f>', if)', s' + , s' x ) denote those values of "a* and 5 which maximize p 2 for some 
source parameters identified by primes to differentiate them from the true source parameters. The following algorithm 
is then defined: 

1. Pick a set of trial source parameters 9' 7 <f)' 7 ip', s' +7 s' x . 

2. Compute ~ct m (6',<f/,i/)',s' + ,s' x ) a,nd~t m (6',(p',il)',s' + ,s' x ). 

3. Form the synthetic response Y from ~ct m and "? m . 



5 



4. Estimate P{Y) using a single detector algorithm. 

5. Retain the source parameters 6', cf>' , ip', s',, s' x if [P(Y) — E[rj\]/E[ri] is the largest to date. 

6. Go back to Step 1. 

The expectation of [P(Y) — E[t]]]/ E[rf\ is just p 2 , so on average this algorithm will converge to the true parameters 
of the source. 

The maximization problem for p 2 can be recast as the maximization of £ subjected to the constraint that E[rf\ is 
constant. The normal equations are 

N 

a j [F t +FfR' ++ (S ij ) + F l +F*R' +x (S ij ) + 

F*F/# x+ (<%) + F*F*R' X x (%) - 
F+F+Ri ++ {5 3l )-F*F+R' +x {5 3l )- 
F+F* R' x + {8 Tl ) F* F/ R' x x (5 Jt )] = (10) 
\ ai a 2 + ai [F+F+s + ■ s + + (F+F* + F*F+) S+ ■ s x + F*F*s x ■ s x ] + 

N 

]T "J [F+F+R++ (% ) + F+F*R +x (% ) + 

F*F+R x + (5 l3 ) + F*F*R xx (5 l3 )} = 0, (11) 

where A is the Lagrange parameter for the constraint, R'^ix) = dRij(x)/dx, and i = 1,2, 3. 
Eq. 1)11) [) can be simplified using the identity Rjj(x) — Rji(-x): 

N 

aj [F+F+R' ++ (S t] ) + F+F*R' +X ($ tf ) + FfF+Xx+fa) + F* F/ R' x x = 0. (12) 

If R+x (x) has an extremum at x = 0, i.e. if R' +x (a;)|a;=o = 0, then a solution to Eq. i|12fl is Si ~ — for i = 1, TV, 
since and R X x(x) are maximal at a; = 0. Back into Eq. (|11J) . this solution gives 

AT 

Aa^ 2 + ^ a,-[F+F/ S+ • s+ + (F+F/ + F^F/)^ • s x + F?F*s x ■ s x ] = 0. (13) 

The choice of if to maximize the signal-to-noise ratio then depends on the angles 9, 0, tp as before, but now only 
on the two numbers |s+|/|s x | &nd s + • s x /|s+||s x | instead of the full waveforms for the two polarizations. These 
numbers are denoted A + / x and A_|_. x , respectively. Eq. I|13|) is an eigenvalue problem (with the Lagrange parameter 
A playing the role of the eigenvalue, and the weight vector ~ct that of the eigenvector), and is straightforward to solve 
numerically. The matrix of the eigenvalue problem is hermitian, so all eigenvalues are real; the eigenvector which 
maximizes Eq. (|5j) with 5y = is picked to form the synthetic response. 

To recapitulate, if i?+ x has an extremum at zero, a synthetic waveform can be constructed from the data of all 
interferometers such that this waveform is optimal for its processing by power detectors, and the maximization of 
the detection statistic for parameter estimation can be performed over only four parameters: 9, <fi, A + / x and A + . x . I 
do not include ip in this list of parameters because it is completely degenerated with A + / x and A + . x . In particular, 
changing the basis where the two polarizations are defined by a rotation (i.e., a redefinition of the polarization angle 
ip) does not change the value of R' +x (x)\ x= q. For the reminder of this paper, I will take tp = 0; this fixes the definition 
of s + and s x with respect to the frame of reference of the detectors network. The CPF algorithm is then defined as 
followed: 

1. Pick a set of trial source parameters 9', 0', A^, x , and A' + . x . 

2. Compute if from Ea. (|13fl . and set (5=0. 

3. Form the synthetic response Y using Eq.©. 

4. Use a single interferometer power detector to calculate P = \Y\ 2 . 
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5. Retain the source parameters 8', <f>' , A^/ x , and A' + . x if P is the largest to date. 

6. Go back to Step 1. 

If the signal is linearly polarized, it can be written as s + = s cos 2tp and s x = — ssm2ip, for some polarization 
angle ip and some waveform s. The eigenmatrix in Eq. I|13(l is then of rank one, and consequently only has a single 
non-trivial eigenvalue. The corresponding eigenvector is given by <n — [F^(6, <fi, 0) cos 2ip — i 7 ^* (#,</>, 0) sin 2-0] /of = 
Fj~(0,(f>,ip)/(Tj for i = 1,2,3. Hence, when the GW signal is linearly polarized, the signals from the different inter- 
ferometers are optimally combined by weighting them with the ratio of their noise variance to the value of the beam 
pattern functions, weighted appropriately by the polarization angle. The signals from the different interferometers of 
the network are therefore emphasized linearly in the observed power of the GW signal, and inversely in the power of 
the noise. 

Another interesting subcase involves directed searches: in that case, the position [6, </>) of a potential GW source is 
known precisely (5^ = 0), and the goal is to be maximally sensitive to gravitational radiation from that source. For 
hi = ai<7i, the maximization of the signal-to-noise ratio can be rewritten as the maximization of the quadratic form 



N 

= J2 MjMij (14) 

<>3=1 



subjected to 



N 



= i, (is) 



where the matrix M has elements M« given by 



» f |S+ll s x| 



(Tier- 



F+F+A +/x + (F+F*+F*F+)A + . X + 

A +/x 



(16) 



It is well known from Rayleigh's Principle [2(j that the maximum value of p 2 is given by the largest eigenvalue of 
the matrix M. This maximum signal-to-noise ratio can be compared to the signal-to-noise ratio that can be obtained 
using the best interferometer in the network, which is typical of the signal-to-noise ratio of an incoherent search. It 
is given by 

Pbcst = max M u . (17) 

1—1...N 

Considering the HLV network simplified so that all interferometers have the same noise level (of ^constant), and 
fixing the beam-pattern functions by selecting a source along the northern hemisphere normal of the HLV plane, the 
ratio p/pbcst varies between 1.03 and 1.57, depending on the value of A + / x and A + . x . Figurc^shows that variation 
of the ratio p/pbest with A + / x and A + . x : when the signal has a large degree of linear polarization, the improvement 
is fairly large, with p/pbest ~ 1.4. When both polarizations roughly contain the same power, the improvement can 
be large 1.5) or very small (~ 1), depending on the structure of the signal (i.e., on A+. x ). If we add a detector at 
TAMA's location to the HLV network, with noise similar to the noise of the other interferometers, the ratio p/pbest 
varies between 1.15 and 1.79 for a source at the same location as before. The four interferometer network works 
better than the three interferometer HLV network, although the improvement is somewhat limited by the unavoidable 
misalignment between instruments located on different continents. 

If R+x does not have an extremum at zero, the normal equations are not necessarily satisfied at Si — — Aj. They 
become non-linear and rather complex to solve, but, more importantly, a full knowledge of the waveforms for the two 
polarizations is then necessary to obtain a solution. Physically, this is a result of the fact that the two polarizations 
may interfere together constructively when shifted by a non-zero lag; this lag gets added to the estimated time delay 
between the interferometers in the network, and a systematic error in the source position estimate appears. It is 
important to realize that this systematic error does not significantly reduce the detection capabilities of CPF, but only 
its positioning ability. Essentially, the convergence to a secondary maximum provides an alternative way to cross the 
detection threshold when the global maximum of the signal power is plagued by noise, such that the measured power 
at its position is not a global maximum of the measured network power. It is nevertheless interesting to estimate the 
size of the position systematic error when the CPF algorithm is used on a signal for which R' +x (x)\ x= o — does not 
hold, i.e. when Eq. (jT3|> is imposed as a solution. This is the subject of the following subsection. 
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The ratio p/pbcst as a function of A + / x and A+. x , for a source along the northern hemisphere normal of the HLV 



Before this question is addressed, however, it is worth examining the physical meaning of the condition 
Rf +X (x)\ x —o = 0. As it was pointed out by [2ll Eqns. 2.30d, 4.3, and 5.18a], the GW strain for slow-motion, 
weak gravity sources radiating mostly due to variations of their mass moments (by contrast to current moments) can 
be written as 



OO I .J - 

F( t )«^^(vv^K/ 

1=2 m=-l J 



p(t - r/c)Y 



lm *r 2+l drdtt, 



(18) 



where STT means "symmetric transverse-traceless" , Y lm are the spherical harmonics, and p is the mass density of 
the source. A number of possible GW sources consist in anisotropic mass distributions that are rapidly rotating 
about a well-defined axis, and consequently radiate principally with I — 2, m = 2 (e.g., binaries, bar or fragmentation 
instabilities, longest live mode of a perturbed Kerr black hole, as noted by j2^). The imaginary part of Y 22 is rotated 
by 7r/4 about the polar axis (the rotation axis) with respect to its imaginary part. This results in the mass distribution 
being sampled at any given time according to two spatial patterns rotated by 7r/4 with respect to each other. Since 
the GW are emitted at twice the rotation frequency of the source, this rotation angle produces a tt/2 phase difference 
between the real and imaginary parts of the time dependent integral in Eq. Ijl8|) ■ As noted by [22j , the dependence 
of the polarization on the inclination angle comes from the pure-spin tensor harmonics (VVl^ m ) STT , which for I = 2 
can be found in |23|. For m = 2, 



(vvr 22 ) STT oc (i 



cos' 



2i cos i e> 



(19) 



where i is the inclination angle of the rotation axis with respect to the line of sight (i — along the polar axis), and 
where e + and e x are the unit linear-polarization tensors for the GW. Consequently, the plus and cross polarization 
waveforms can be written as 



h+ cx (1 + cos 2 i) cos <&(£) 
h x oc 2cosisin$(t), 



where $(t) is some phase function. 

By definition, equations (|2C)|) and (|21|l give 



/OC P oc 
/ dtidt 2 Q{ti,t 2 ) cos <$>(ti) sin <&(t 2 - r), 
-oo J — oo 



(20) 
(21) 



(22) 
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so that 

/OO POO 
/ dtidt 2 Q(ti,t 2 )<f>' (t 2 ) cos $(tx) cos <Z>(t 2 ). (23) 
-oo ^ — oo 

In general, the phase will have at least a linear component ~ wi], so that 

R' +X (x) \ x= q oc — cos l(1 + cos 2 l)uj. (24) 

As the system's polar axis becomes aligned with the line-of-sight, the assumption that R' +x (x)\ x= o = becomes 
progressively worst, and significant systematic position errors might appear, and will have to be accounted for, as 
discussed in the next section. Some signals might also satisfy the condition R' +x {x)\ x =q = 0. This is the case, for 
instance, if Q(t) is an even function of t. 

Sources which radiate predominantly in a I = 2, m — 1 mode will show a similar correlation between their plus and 
cross polarizations. In that case, the real and imaginary parts of Y 21 differ by a ir/2 rotation about the polar axis, 
but have a m = 1 symmetry, so that the waves are radiated at the spin frequency, and consequently the phase shift 
between the plus and cross polarizations is again n/2. It might be that the dominating population of sources to be 
observed will not be dominated by rotation about a principal axis; as a result, the correlation between the plus and 
cross polarizations might be quite arbitrary. In an axisymmetric core collapse, for instance, the I = 2, m = mode 
dominates [24|. The angular response is (VVY 20 ) STT oc sin 2 t e + , so the waves are linearly polarized. 



A. Systematic Position Errors 

Eq. I|12fl can be used to check the error on 9 and </> by solving it for 5 , with it obtained from the CPF algorithm, i.e., 
from the solution of Eq. 1)13(1 . For every trial choice of (6 1 , 0', A' + ^ x , A' + . x ), the algorithm returns a set of weights it 

(by definition of CPF, 5 = 0). However, the maximum of the signal-to- noise ratio occurs when the normal equations 
are satisfied; assuming that we let A + / x and A + . x be varied freely in the maximization of P, the values of 9 and <f> 
returned by this maximization will be those that solve Eq. (|12ll . 

Consider the following parameterization of the signal cross-correlation functions: 

R' ++ {t) = R' xx {t) = -uj 2 t + 0^) (25) 
R' +x (t) =uj a - u\t 2 + 0(t 4 ) (26) 
R' x + (t) = -coo + u\t 2 + 0(t 4 ), (27) 

for some parameters loq, oj\, and lo 2 - If s+ and s x were long, nearly monochromatic signals of angular frequency 
u>, with the same amplitude but a phase difference of tt/2, for instance, the parameters could be chosen to be 
luq = lui = lo 2 = lo. In order to get an analytical solution, I linearize Eq. I|12l) to obtain the first order equation 

E a j [F ij u, 2 1 sl?-u; G ij }=0, (28) 

where 5y represent first order errors between the true time delays and the delays returned by the algorithm, and 
where = FfF* + F* F+ and G i:j = FfF* - F*F+. The solutions to this linear system of equations are 
degenerate; for N = 3, they are 

^i ^(i) _ a<iFi 2 Gi3 + a 2 Fi 2 G 2 3 + a 2 F 2 ?,Gi 2 + a^F^Gi^ 
coo 13 aiFi 2 Fi3 + a 2 Fi 2 F 2 z + 03^13^23 

and 

&x ^(l) _ aiFuGis + a 2 Fi 2 G 2 3 — aiFi^Gi 2 + a3Fi3G*23 
coo 23 ai-Fi2-Fi3 + a 2 Fi 2 F 2 3 + a^Fi^F^ 

The times-of-flight between two pairs of detectors are sufficient to triangulate the position of the source on the 
sky and to obtain its position (8, </>), up to a reflection with respect to the plane containing the three detectors. The 
magnitude of the systematic error due to R' +x (x)\ x= o 7^ is given by the arclength of the portion of the great 

circle connecting the true source position and the position obtained by triangulation from and in Eqs. I|29|) 
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and (|30|) . To give a representative example, is computed for the long monochromatic signal described above, for 
the HLV network with af — 1, i.e. under the simplifying assumption of instruments with identical noises at all sites. 
Only a very small fraction of the sky (about 0.7%) has a negligible systematic error [l^(2n x 40 Hz/a;) < 0.01 rad]. It 
is therefore plain that position estimations will be grossly off target if the assumption that R' +x (x)\x=o = is wrongly 
made. 

When R' +x (x)\x=o ^ 0, more information about the waveforms is required to estimate the source position. One 
possible approach is to use the CPF algorithm, which assumes R' +x (a;)|x=o = 0, and then to correct the source position 
estimate using Eqs. (I29[) and (|30|1 . This requires only a knowledge of the slope of R+ x with respect to the slope of 
R++ at zero lag (i.e., ujI/ujo), which is closely related to the characteristic frequency of the waveforms. A map from 
estimated position (with systematic error) to actual position must be constructed for every choice of lo\/ujq. The 

(2) 

remaining systematic error is given by the higher order terms not included in the correction. For . the second order 
errors on the time delays, the second order equation derived from the linearization of Eq. (|12[1 is 

N 



E 



- "iGiMf + 2<$ ) <5g ) )] = 0. (31) 



This linear system of equations can be solved to obtain: 

^ifr^ = a 1 S^[2a 2 s[ 2 ) F 12 G 12 + a 3 S^G 13 (F 12 - 2^S^G U )] + a 3 {a 3 S$ 5$ G 13 [F 23 + 2^sg 3 ] G 23 ] 
6\ 3 > ^1 w i 



a 2 [Sg )2 F 12 G 23 + 5$G 12 (5$F 23 + 2^sg\dg ) + 5$)G 23 )}} (32) 



and 



A 2 ) r> 33 

-4fr~ = a 3 Sg 3 ] G 23 (F 12 + 2^d[ 1 2 ) G 12 ) + oa^fFia + 2^ 1 3 ) Gi 3 )] + ax{2a 2 8^ 8^ F l2 G 12 
6 23 > w i w i 



where 



uj 1 



D = a 3 ^{ ai sg\F 12 - 2^5$G 12 )(F 13 + 2^ 1 3 ) Gi 3 ) + 

0^2 ^1 ^1 

0F 23 + 2 t 4sV 3 ) G 23 )[a 2 5{ 3 ) (F 12 + 24^2^12) + a 3 s{\\F 13 + 2^6^ G 13 )]}. (34) 

LiJ-y UJ^ UJ-I 

For the arclength of the portion of the great circle connecting the true source position and the position obtained 
by triangulation using i.e., for the systematic error after correction using the first order expansion, and for 
the same example as above (with u> = 807r rad/s), one finds that 23.2% of the sky has a negligible systematic error 
(l^ < 0.01 rad). Figure shows the fraction of the sky with a systematic error smaller than a certain value, for 
and /( 2 ). It should be noted that for only about 50% of the sky is it meaningful to use an expansion in the small 
parameters 6^ and 5y . All in all, these numbers show that it is possible to use the CPF algorithm for any signals, 
at the cost of limiting the sky coverage to ~ 25%, and of requiring one additional piece of information about the 
signal, the value of the ratio uj 2 /ujq, which describes the behavior of the signal cross-correlation functions near zero 
lag. Shorter signals, or signals with less overlap between their two polarizations, are likely to offer smaller systematic 
errors on position estimates, so that the correction suggested above becomes unnecessary. 



IV. NUMERICAL SIMULATIONS 



All the numerical simulations were performed using the CPF implementation for the LIGO Data Analysis System 
(LDAS), which uses the TFCLUSTERS algorithm 0] to measure the signal power in a given time-series. The HLV 
network was analyzed, except that the noises of all three interferometers were assumed identical, for simplicity. For 
every realization of the simulation procedure, a 10 s long segment of white noise of unit variance, sampled at 16384 
Hz, was generated independently for all three interferometers. The use of white instead of colored noise is not an 
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FIG. 2: The fraction of the sky with a systematic error smaller than the value of the position error plotted on the horizontal 
axis, for a long monochromatic signal of angular frequency lo = 807r rad/s, with two polarizations of the same amplitude, but 
at a phase offset of 7r/2. The rightmost curve is for the leftmost one for l^ 2 \ 

important loss of generality, since TFCLUSTERS is quite robust with respect to the presence of correlations in the 
background noise. 

Simultaneously, two 1/8 s long segments of unit variance white noise were generated, and filtered by a 6 th order 
elliptical bandpass filter with 3 dB cutoff frequencies at 125 Hz and 150 Hz, so that the amount of power outside this 
band was negligible after filtering. These two segments were then truncated to their central 1/16 s long portion, and 
were used as the plus and cross polarizations of the gravitational wave signal incident on the network of detectors. By 
construction, the average value of A + / x was 1, and that of A + . x was zero. In some simulations, the plus polarization 
waveform was multiplied by \pl and the cross polarization waveform was divided by y/2, so that on average A + y x 
was equal to 2. 

All gravitational wave signals were injected from the position corresponding to the northern hemisphere normal of 
the HLV plane, so that they arrived in phase at the three interferometers. The two polarizations were combined at 
each interferometer using the beam-pattern functions and F* , and were added directly to the background white 
noise. 

A scale factor A was multiplying the GW signal injected in the simulated data for all interferometers. If the 
frequency band occupied by the signal, its duration, and its arrival time had all been known exactly, it would have 
been possible to filter the data optimally in that band. In that case, the signal-to-noise ratio p pt for a certain scale 
factor A would have been given by 



X^ + O- (35) 
i=l 

Of course, it cannot be assumed that this information is available, but I will nevertheless use Eq. l|35|) as a measure 
of the strength of the injected signals, instead of using A which is less intuitive. Numerically, p op t(A) = 1.34A for a 
source at the northern hemisphere normal of the HLV plane. 

The simulations were performed using a realistic hierarchical implementation of CPF. The simulated data from the 
three interferometers were first processed separately by TFCLUSTERS to produce three lists of events. In terms of the 
notation developed in the settings of TFCLUSTERS were a = 0, a = 5, and S = [0, 0, 0, 0, 0, 0, 2, 3, 4, 4]. Only the 
frequencies below 1024 Hz were considered, and the time resolution of the time- frequency decomposition was T =1/8 
s. The number of events in this first stage was controlled by the black pixel probability threshold, po (note that larger 
values of po give larger false alarm rates). 

The events produced by TFCLUSTERS are in the form of rectangles in the time-frequency plane, with information 
about the power present in each pixel of these rectangles. For a given cluster identified by TFCLUSTERS, the rectangle 




p pt(^) = A 



\ 
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is defined as the smallest rectangle containing all the pixels of the cluster. A triple-coincidence condition was therefore 
applied as in |25|: to be coincident, the time-frequency rectangles corresponding to the events from all three inter- 
ferometers had to be overlapping. This coincidence condition selects events that are close in time and in frequency, 
and can be understood as a standard time and frequency coincidence gate with varying windows that are fixed by the 
events under consideration. 

All the coincident events that were present in a given 10 s segment were then considered in turn. Their start time, 
duration, central frequency and bandwidth were estimated from the smallest rectangle in the time-frequency plane 
that could contain the union of the rectangles from the individual events. The CPF algorithm was then ran on the 
data, once for each coincident event. The implementation employed used TFCLUSTERS to process the synthetic data, 
and the measure of the power was the sum of the power in all the pixels identified by TFCLUSTERS that were inside 
the time-frequency rectangle identified from the triple coincidence. The parameters for TFCLUSTERS were the same as 
those mentioned above, except for the black pixel probability, which was set to a value p%. The threshold defined by 
Pi was chosen so that only loud enough signals were detected by TFCLUSTERS, and their estimated power was linearly 
related to their actual power. 

The power measured by TFCLUSTERS was maximized over the source position (two angles) and over the parameter 
A + . x . It was assumed that the value of A + / x was known beforehand, in order to keep the size of the parameter space 
small enough for simulations. In a first time, the sky was covered by picking 100 points uniformly distributed in the 
range [0, 2n[ for the right ascension, and 100 points in [—1, 1[, uniform in the sine of the declination angle. In addition, 
10 points were used to cover the range [— 1, 1[ uniformly for the parameter A + . x . Consequently, TFCLUSTERS was ran 
10 5 times on every 10 s long simulation. Including the overhead from the LDAS system, this part of the search ran 
in ~ 225 s on 31, 2 GHz Pentium IV computers, with 512 Mb of RAM. 

If none of the triple coincidence events registered above the pi threshold of TFCLUSTERS when analyzed by CPF, a 
non-detection was reported and the analysis was stopped. Otherwise, a detection was announced, and CPF produced 
a scan of the parameter space for every triple coincidence event above threshold. The triple coincidence event with 
the largest maximum power was then selected as a possible GW candidate, and was analyzed in more details. The 
point in parameter space where the power was maximum defined the parameters for a second run of CPF, used to 
obtain a refined position estimate not limited by the coarseness of the grid covering the parameter space. The value 
of A+ . x was fixed to the value estimated in the first run, and a square search window of size 0.2 rad in right ascension 
and in declination, with 50 steps in both angles, was centered on the value of the position obtained in the first run. 
The position with the maximum power in this second run was taken as the final estimate of the source position. 

A. Detection 

For the sole purpose of detecting the presence of a signal in the data, only the first run of CPF was required. Data 
were ran through TFCLUSTERS separately, the events were combined in the triple coincidence gate, and the data were 
fetch through CPF for the parameters defined by each triple coincidence event. If at least one of the triple coincidence 
event lead to a detection by CPF, the 10 s long time interval under scrutiny was assumed to contain a signal. By 
design, this hierarchical scheme required a fairly permissive threshold in the first stage where TFCLUSTERS was ran 
independently on every interferometer, so that a given signal was very likely to make it to the second stage where 
CPF was ran. The limit on this threshold was determined by the availability of computational resources, and by 
the confusion that resulted from the proliferation of events at low threshold. Most of the rejection of accidental 
coincidences occurred at the second stage, where CPF was operated with a reasonably strict threshold. 

A numerical experiment was performed by running this simulation a large number of times, with and without signal 
injection. When injected, the signal had p opt = 13.4 and A + / x = 1. The thresholds were chosen to be po = 0.14 
and pi = 0.012. In the first stage (triple coincidence), the probability to detect a signal was Pd = 0.92 ± 0.01, and 
the probability of a false alarm when no signal was present was Pp — 0.62 ± 0.02. In the second stage (CPF), it 
was measured that (Pd,Pf) — (0.86 ± 0.01, 0.10 ± 0.01). Overall, it was measured with both staged combined that 
(Pd, Pp) = (0.79±0.02, 0.064±0.009). Table|I|gives the detailed results from the simulations. The errors quoted here 
come from 68.3% confidence intervals ("lc") for a Bernoulli process, built using the prescription of [2(|. The 6.4% 
probability of false alarm is larger than typical values for GW searches, as it would give a false alarm rate around 7 
mHz. It was chosen, however, in order to provide enough detections for small errors on the measured probabilities. In 
a more realistic setting, po would be similar to the value used here, while p\ might be smaller by one or two orders of 
magnitude. It should be noted that the threshold settings were found by trial and error; there is an infinity of points 
along a curve in the pojPi _ pl a ne that give a 6.4% probability of false alarm, and a choice different than the one above 
may give a larger probability of detection. 

Nevertheless, it seems that the choice I made for the po and p\ thresholds is sufficient to show the superiority of 
the coherent approach over the incoherent one for detection. I repeated the experiment above, but using only the 
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no signal 


injected signal 




injected 


with p pt = 13.4 


total number 


981 


730 


detected by 


611 


669 


triple coincidence 


[0.606,0.640] 


[0.904,0.927] 


undetected by triple 


370 


61 


coincidence 


[0.361,0.394] 


[0.0726,0.0962] 


detected by CPF 


62 


578 




[0.0549,0.0731] 


[0.775,0.808] 


undetected by CPF 


549 


91 




[0.542,0.576] 


[0.112,0.139] 



TABLE I: Details of the simulations to measure the efficacy of the hierarchical implementation of the CPF. Numbers in brackets 
show the 68.3% ("If") confidence interval for the fraction of the number of trials to the total number of trials. 



first stage to detect events. A detection was announced when at least one triple coincidence was observed between 
the outputs of the three TFCLUSTERS runs on the interferometers' data. In that case, a threshold of po = 0.11 gave 
(Prj, Pf) = (0.70 ± 0.02, 0.07 ± 0.01), in an experiment with 742 trials for the measurement of Pd (519 detections), 
and 742 trials for Pp (48 detections) . At a similar false alarm probability, the probability of detection is significantly 
smaller in the incoherent case than in the coherent case. For the particular signal and false alarm probability under 
consideration, the signal-to-noise ratio would have to be increased to p op t = 16.8 in order for the incoherent approach 
to be as efficient as the coherent one. For a homogeneous distribution of sources in space, this corresponds to a factor 
of <~ 2 improvement in detection rate, assuming no significant degradation of the CPF algorithm performances with 
respect to the incoherent algorithm as the position of the source is varied away from the northern hemisphere normal 
to the HLV plane. 

The performances of the CPF search were mostly limited by the quality of the estimation of the time-frequency 
rectangle containing the burst, in the first stage of the analysis (the triple coincidence). With the signal injection 
for p opt = 13.4, the first stage gave an estimated rectangle that overlapped with the one containing the signal in 
96% of the 578 cases where the signal was discovered, but only in 3% of the 61 missed detections was this the case. 
Without signal injection, only 1% of the 611 triple coincidences had time-frequency rectangles overlapping with the 
signal rectangle. These numbers show that if an oracle were available to provide the rectangle containing the signal 
without error every time the search is ran, a probability of detection > 97% would be possible for a probability of false 
alarm < 1%. Stated differently, CPF is extremely efficient at detecting a burst when it receives the right parameters 
describing that burst; the triple coincidence incoherent search provides many candidates; when one such candidate 
corresponds to the signal, CPF picks it out of the others very efficiently. A better approximation to this oracle than 
the one used here might be to tile the time- frequency plane with a variety of rectangles, and to run CPF on each 
rectangle. 

Suppose that it is known that the signal has a bandwidth of 25 Hz. One can cover the time- frequency plane with 
non-overlapping rectangles of duration 1/16 s, and bandwidth of 25 Hz, so that the 10 s long, 1024 Hz bandwidth 
data segment in my simulations is covered by ~ 6500 tiles. The false alarm probability for each run of CPF must be 
reduced to ~ 10~ 5 so that the global search has Pp ~ 7%. The efficiency of the search would then approximately be 
given by the probability of detection of CPF with an oracle, for a threshold giving Pp = 10~ 5 . I have not measured 
this probability of detection, but it is plausible that it is larger than the 79% efficiency measured for the hierarchical 
implementation used in the simulations. However, running this search in real time would be prohibitively expensive, 
as it would require ~ 10 4 teraflops of computational power. In the present hierarchical implementation, the first stage 
took <~ 30 gigaflops to run in real time, and the second stage, ~ 800(Pf/0.62) gigaflops, where Pp is the false alarm 
probability in the first stage. These numbers should be taken as upper bounds on the required computational power, 
because the codes were not optimized to minimize overhead, to make an optimal usage of the parallel resources, or 
to do an optimal scan of the parameter space. Optimized codes should be able to run at least <~ 5 times faster. 

B. Position estimation 

In order to get a precise idea of the magnitude of the position estimation errors, the analysis scheme described above 
was simplified by removing the first incoherent step involving the three different instances of TFCLUSTERS running at 
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FIG. 3: Scatter plot of the estimated position of a source injected along the northern hemisphere normal of the HLV plane, 
from 240 realizations of the simulation with A + / x = 1 and p op t ~ 13.4. Horizontal axis is right ascension, vertical axis is the 
sine of the declination. The curves represent loci of equal time delay for the three independent interferometer pairs. The point 
where they intersect in the upper left corner of the figure is where the signal was injected. 

each site. Instead, the CPF algorithm was instructed to compute the power according to the output of TFCLUSTERS 
in a rectangle of duration 1/8 s located at the right position in the time series, with a lower frequency of 50 Hz and 
an upper frequency of 150 Hz. The black pixel probability was set to p\ = 5 x 10 -3 , so that the number of clusters 
unrelated to the signal and produced only by the noise was small. 

Figure [3] presents a scatter plot of the position estimates obtained for 240 realizations of the simulation, when 
A + / x = 1. The estimates tend to cluster along the curve corresponding to the locus of positions having equal delays 
at the Hanford and Livingston interferometers. This is a direct consequence of the good alignment between these two 
interferometers, and the relatively poor alignment of the Virgo detector with them. Estimates tend to fall on that 
curve, but also to cluster at different places along it, where the signals at Hanford and at Livingston are delayed by 
an integer number of the characteristic periods of the signal with respect to the signal at Virgo. 

Figure 0] shows a similar plot as FigOU except that A + / x = 2. This corresponds to a GW signal which has more 
structure in its polarizations than the one for the case A + / x = 1, i.e. which is closer to a linearly polarized signal. 
Since linearly polarized signals are the easiest ones to analyze with a network of interferometers, it is expected that 
the position estimates will be better. This is indeed was is observed in FigQ] the estimates still hug the Hanford- 
Livingston equal-delay curve, but now present less scatter around the points where the signals are in phase with the 
Virgo signal along that curve. 

The position error can be quantified as in Fig|21 it is taken to be the length of the shortest portion of the great 
circle joining the estimated position and the true source position or its mirror image with respect to the HLV plane. 
The error can be defined with respect to the mirror image because there is a natural ambivalence in the estimation of 
the position when only three detectors are used. Figure [5] shows the cumulative distribution of the position error from 
the simulations, for A + / x = 1 and for A + / x = 2, and for two different values of the signal-to-noise ratio of the injected 
signals. As expected, signals with larger values of the signal-to-noise ratio or of A + / x lead to smaller position errors. 
The curves in Fig[S] present a number of "steps", which are produced by the clustering along the Hanford-Livingston 
equal-delay curve at positions in phase with Virgo. 

Roughly 50% of the trials lead to unusable position estimates (errors > 10 degrees) when A + / x = 2, while this 
number reaches ~ 80% when A + / x = 1. However, with A + / x = 2, approximately 25% of the trials have errors smaller 
than one degree. Moreover, at least in the regime of signal-to-noise ratios under consideration here, the scaling of the 
position error with the signal-to-noise ratio is rather weak. 
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FIG. 4: Same as Figure El but with A +/x = 2. 
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FIG. 5: The fraction of all simulations that gave a position error smaller than the value plotted on the horizontal axis. The 
two continuous lines correspond to A + / x = 1, and the two dotted lines to A + / x = 2. In both cases, the rightmost curve is for 
Popt = 13.4, and the leftmost one is for p op t = 35.6. Each curve is built from 240 realizations of the simulation. The error on 
the curves is estimated to be ~ 5% of the fraction of trials. 

V. CONCLUSION 

A method was presented for the optimal generalization of the power detectors developed for single interferometers 
so that they can process coherently data from a network of interferometers. The coherent method, as compared to 
an incoherent approach where event lists independently generated at all interferometers are searched for coincidences, 
offers the advantage of better detection efficiency, and the possibility to accurately estimate the position of the source. 
A few systematic effects affect the performances of the CPF algorithm for position estimation, including: cross terms 
between the plus and cross polarization waveforms of the GW signal, lack of differences between the characteristics 
of the two polarization waveforms, and misalignment of the interferometers of the network. 
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The three effects are related, and reflect the obvious fact that GW signals incident on a network of misaligned 
detectors will only show entangled versions of their two polarization waveforms, different in each interferometer. If 
one of the two polarizations is significantly stronger than the other (i.e., the GW signal has a stronger degree of linear 
polarization), the problem is drastically simplified. Similarly, aligned interferometers are much less sensitive to this 
problem. If the two polarization waveforms are fairly coherent with each other and of similar amplitudes, the cross 
terms between the two polarizations in different interferometers may show significant maxima when the time shifts 
imposed on the different data streams do not correspond to the differences in time-of-arrival of the GW signals at 
each interferometer. These maxima cannot be distinguished from the maximum resulting from the product of the 
waveforms of the same polarization in different interferometers, and systematic position errors may result. It should 
be noted, however, that the same effect leads to similar position errors when the differences in arrival time of the GW 
signal in different interferometers are estimated by maximizing the cross-correlation between pairs of interferometers, 
and these time differences are used to triangulate the source position. 

The study of the CPF algorithm presented here was based on a three interferometer network. However, the structure 
of the algorithm allows for any number of interferometers. More complex networks will most likely reduce the effects 
of the systematic errors, by increasing the number of linear combinations of the two polarization waveforms that are 
being sampled, or, equivalently, by increasing the area of the sky where at least three interferometers are measuring 
similar combinations of the two polarizations. Numerical simulations on the network formed by the LIGO Hanford, 
LIGO Livingston, and Virgo interferometers (the HLV network) with a short random signal showed that the CPF 
algorithm could be used to accurately measure the position of the source a significant fraction (~ 25%) of the time, 
for reasonably strong sources located at the normal of the HLV plane. In about ~ 90% of the trials, the CPF algorithm 
correctly placed the source position at a point where the signals at the Hanford detector and at the Livingston detector 
were in phase. However, due to its misalignment with respect to the LIGO detectors, the information provided by 
the Virgo detector was often insufficient to pinpoint correctly the position of the source. The average position error 
was a rather weak function of the strength of the signal, but a stronger function of the amount of difference in the 
structure of the two polarization waveforms. 

It was also shown that the CPF algorithm offers better detection efficiencies than its incoherent equivalent, both for 
directed and for all-sky blind searches, and independently of the systematic errors affecting the position estimations. 
In the former case, improvements in the detection signal-to-noise ratio of 40% or better are expected, excepted for 
a few values of the source parameters. In the latter case, a 25% improvement in signal-to-noise ratio was measured 
for typical source parameters. This improvement for the all-sky blind search comes at the cost of increasing the 
computational power required to perform the data analysis in real time by a factor of ~ 5 — 30. This may be 
significant, especially since lengthy simulations probably have to be completed in order to estimate the false alarm 
and detection statistics of a real search. However, the signals studied in the numerical simulations, which could be 
viewed as very rough approximations to the signals that could be emitted in the collapse of the core of a star in 
supernova explosions, would have been detected by the CPF algorithm at a rate ~ 2 times larger than the rate for its 
incoherent equivalent, assuming a homogeneous spatial distribution of the sources, and that the performances of the 
CPF algorithm for a source injected at the normal of the HLV plane are characteristic of its performances at other 
injection positions. 

Finally, an important advantage of the design of the CPF algorithm is that it inherits the robustness, efficacy, and 
computational efficiency of the single interferometer power detectors. In other words, the CPF algorithm should not 
be more sensitive to non-Gaussian or non-stationary noises than incoherent analyzes, it should be searching more 
efficiently for GW signals sharing similar morphologies, and, while more computationally intensive, its speed should 
scales the same way with improvements in the implementation of the single interferometer power detectors, or with 
computing hardware ameliorations. 
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